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We study the effects of gluon interference in the production and semi-leptonic 
decay of a tt pair above threshold at the Next Linear Collider (NLC). We calculate 
all matrix elements to next-to-leading order and use the resulting expressions for 
the development of a Monte Carlo event generator. Our results show effects at the 
level of 10% in differential cross-sections. We thus extend previous results obtained 
by analytical means in the soft-gluon limit. 



Due to its large mass, and, consequently, its large decay width (m t = 176 
GeV, r t = 1.57 GeV), the top quark decays before it has a chance to hadronize. 
It is therefore possible to perform a precision analysis of top quark production 
and decay, which facilitates the identification of effects attributable to new 
physics beyond the Standard Model. The best laboratory for such explorations 
is an e + e~ collider, because the initial state is free of color effects. Such 
a machine (the Next Linear Collider (NLC)) will hopefully be built in the 
not-so-distant future. Detailed studies of processes through which one may 
extract useful information have already been performed El. Gluon radiation 
contributes at the 10% level, in general. Although it does not modify inclusive 
cross-sections, it has a non-negligible effect on various final-state distributions. 
Therefore, it has to be included in a precision study of top quark physics. In 
studying gluon radiation, one cannot ignore the effects of gluon interference. 
The latter depends upon the value of T t and may be an important tool in the 
study of this parameter, in addition to precision studies H. 

We have studied gluon interference effects in the process 

e+e- -> 7* , Z* -> it -> W + b W~b g -> i + vb fvb g (1) 

We have calculated all amplitudes analytically, with the aid of the symbolic 
manipulation package FORMEl. The resulting expressions were then used for 
the development of a Monte Carlo event generator. The generator was based 
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on the C++ program written by Schmidt for the study of the same processta, 
in the narxpw top width approximation. Details of the calculation are reported 
elsewherelx The results from our simulations are in agreement with analytical 
results obtained by Khoze, et al. U in the soft gluon regime. They are also 
similar to the results obtained by Petex. and SuminoEj in the threshold region, 
where color Coulomb effects dominate □. 

At tree level, the amplitude for tt production and decay factorizes. With 
polarized e + e~ beams, one can perform a clean helicity analysis, because the 
spins of the top and anti-top quarks get transferred to the final products. It 
is then advantageous to express the amplitudes in terms of helicity statesQ. 

Loop diagrams at o(a s ) include a pentagon diagram (gluon exchange be- 
tween the b and b quarks), which renders the calculation of o(a s ) effects cum- 
bersome. To systematically calculate these effects, observe that all loop di- 
agrams contain one or more internal t or t legs, which are rapidly decaying 
quarks. To take advantage of this fact, we define the quantities 

rn = (pi - m 2 t )/m 2 t , 772 = (p 2 - m 2 t )/m 2 t . (2) 

The momenta of t and t are defined in terms of the quarks' final decay products. 
rji and 772 are small parameters, o(T t /mt). Physically, this means that t and 
t are produced nearly on shell. It is advantageous to expand the amplitude in 
these parameters, 

+ o« 2 ). (3) 

where A needs to include both virtual and real gluon contributions (soft and 
collinear), in order to avoid singularities. 

After calculating amplitudes analytically at first order in 771,2, we used the 
formulas to construct a Monte Carlo event generator. The generator is based 
on a program written by Schmidt in C++. Schmidt's program generates events 
at tree level without gluon production with 100 % efficiency. This is achieved 
by making use of helicity states, not only for the final states, but also for the 
intermediate t and t quark states. The amplitude is then expressed in terms 
of helicity angles which are defined in different Lorentz frames for the three 
vertices, respectively. This construction is carried over to one-loop level, if one 
neglects gluon interference effects, because the amplitude is still factorizable at 
that level. The virtual, soft and collinear corrections amount to corrections in 
the form factors for the three vertices in the diagrams. Thus, one can construct 
a very efficient event generator. 

To account for gluon interference effects, one may not rely on helicity 
angles, because the amplitudes are no longer factorizable, and need to be ex- 
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Figure 1: Distribution of gluon energy at center-of-mass energy of 400 GeV. In (a), we 
exhibit interference effects. In (b), we show the dependence on IV 



pressed in terms of a single Lorentz frame. This complicates their calculation. 
To simplify the expressions, we have introduced polarization vectors for the 
intermediate boson states defined with respect to the b (b) quarks (which 
are treated as massless). This quadruples the number of diagrams, but sim- 
plifies the final expressions. At tree level, we have checked the accuracy of our 
expressions numerically by using Schmidt's program. At loop level, we have 
performed a similar check in the limit 771,2 — ► 0. 

Since gluon interference effects are relatively small, we took advantage 
of Schmidt's efficient code, by keeping the event generating procedure and 
introducing the new effects as a correction to the weight. 

We set the mass of the b quark to zero in matrix elements. This is essential 
for the simplification of the expressions for the various amplitudes. We also 
neglected effects of o(r]f 2 ), which include gluon interference between the b and 
b quarks. Such effects are very small and are not likely to contribute at the 
desired accuracy for the NLC. Should the need for higher accuracy arise, our 
method can be readily applied to accommodate these effects. We also did 
not couple the program with a jet algorithm. This would be necessary for 
a direct comparison with experimental data. Initial-state radiation effects E3 
were included in the program, but will not be shown here. 

Next, we present numerical results obtained with the Monte Carlo event 
generator. As has been shown, there is a complete cancellation^between vir- 
tual and real gluon contributions to an inclusive cross-section 0. Moreover, 
differential cross-sections receive contributions from gluon interference effects 
to o(a s T t /mt). The major contribution to such corrections comes from soft 
gluons of energy E g ~ T t u. Our results include hard gluon effects and are,in 
general agreement with the semi-classical analysis of soft gluon productionBEl. 
Indicatively, we plot the distribution of the gluon energy (Fig. ||) at center-of- 
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mass energy y/s = 400 GeV. We have introduced a jet resolution parameter 
Vcut — 0.05 by demanding (k + pb) 2 /s,(k + pi) 2 / s > r)^ ut where fc M is the 
four-momentum of the gluon. We have imposed no other cut on the gluon. 

In Fig. ^( a), we have compared distributions with and without interference. 
In general, interference effects enter at a 10 % level, so they need to be included 
in a precision study of top quark production and decay. Moreover, in Fig. Q(b), 
we show the dependence of the distribution on the decay width of the top 
quark, r t . The decay n width has a most pronounced effect in the range of 
gluon energies E g ~ T t □. This shows that gluon interference effects in the soft 
gluon regime can be used to measure the_top quark decay width. This can 
already be seen by a semi-classical analysisa, but the full quantum mechanical 
calculation, which is done here, is needed fat an accurate theoretical prediction. 
Further discussion can be found elsewhereH. 
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